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in the plane of sky) and the observer (assumed to be at z = oo), 
the effect is a plain change of the original polarisation angle xo 
according to 
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ABSTRACT 

We develop a novel maximum-a-posteriori method to measure magnetic power spectra from Faraday rotation data and implement it 
in the REALMF code. A sophisticated model for the magnetic autocorrelation in real space permits us to alleviate previously required 
simplifying assumptions in the processing. We also introduce a way to treat the divergence relation of the magnetic field with a 
multiplicative factor in Fourier space, with which we can model the magnetic autocorrelation as a spherically symmetric function. 
Applied to the dataset of Hydra A north, we find a power law power spectrum on spatial scales between 0.3 kpc and 8 kpc, with 
no visible turnover at large scales within this range and a spectral index consistent with a Kolmogorov-like power law regime. The 
magnetic field strength profile seems to follow the electron density profile with an index a = 1. A variation of a from 0.5 to 1.5 would 
lead to a spectral index between 1.55 and 2.05. The extrapolated magnetic field strength in the cluster centre highly depends on the 
assumed projection angle of the jet. For an angle of 45° we derive extrapolated 36 fiG in the centre and directly probed 16/jG at 50 
kpc radius. 

Key words, galaxy clusters - magnetic fields - Faraday rotation 

1. Introduction 

1.1. Motivation 

Intergalactic magnetic fields are of great scientific interest. They 
are part of the intergalactic metabolism, they affect the propaga- 
tion of heat and cosmic rays, and they are tracers of the dynam- 
ical state of the intergalactic medium. One way to retrieve their 
properties is to measure the Faraday rotation effect they cause 
on traversing polarised light. This has lead to significant insight 
into the presence, strength, and morphology of magnetic fields 
in galaxy clusters. 

Although it is not possible to retrieve the exact magnetic 
field in position space from Faraday rotation data, the magnetic 
power spectrum is accessible and contains information on the 
typical field strength, correlation length, and the distribution of 
magnetic energy onto different length-scales. All these quanti- 
ties can be predicted in magneto-genesis theories. The latter, the 
spectral energy distribution, is a sort of fingerprint of these theo- 
ries. Magneto-genesis-scenarios in which the energy is injected 
on large scales by gas stirri ng typically exhibit K olmogorov- 
like p ower law spectra (Su bramanian et al., 2006; E nBlin & Vogtl 
1200(^1 . Whereas scenarios in which the magnetic fields are built 
up fror n small to large scales typicall y exhibit different spectral 
slopes (ISchekochihin & Cow lev 2006). Magnetic field measure- 
ments by Faraday rotation therefore have the potential to distin- 
guish between these theories. 

1.2. Faraday rotation 

The plane of polarisation of light traveling to us through a mag- 
netised plasma is changed by the wavelength-dependent Faraday 
rotation effect. If the Faraday rotating plasma is entirely located 
between the radio source (at physical depth and at position 



(1) 



with a proportionality constant gq = /(Inmlc*), with c, e, and 
nig being the speed of light, the electron charge, and the mass, re- 
spectively. Furthermore, «e(x) and B(x) are the electron density 
and the magnetic field within the Faraday rotating foreground 
medium, and 0(xx) denotes the Faraday depth of the source at 
x^. Although the original polarisation angle ;ifo is unknown, mea- 
surements of the polarisation angle at a range of wavelengths A 
provide this observationally, and also an estimator of the Faraday 
depth (p via a fit to Eq. [1] The observational estimate of the 
Faraday depth is called the rotation measure (RM), because it 
expresses the rate at which the linear polarisation rotates as a 
function of A^. 

Some of the many st udies of magnetic fields via Farada y 
rotation are iDennisoiil ( 1979h, Lawler & Dennison (1982), 
DreheretalJ (Il987h iLaind ^1988h. iHennessv et al.1 (1 19891) 
Tavlor et al.' (199^), 'Goldshmidt & Reohaeli' (1993'), 'Kim et aL 
1991), Ge & Owen ( 1994), Feretti et al. ( 1995), Johnson et al 

19951). iFeretti et alJ (Il999allbir 
oretalJ ( I2001'). 'Carilli & Tavlor 



i 

(f 19951). iJohnson et all i 
ICIarkeetalJ d2001h iTav 



. [Tavlor et al. 



j20m 



Eilek & Owen (2002), Murgia et al 



(2004), Laing et al. (2006'),'Xu et al.'(2006), Tavlor et al. (200' 
Cho & Ryu '(2009), Feain et al. (2009), and Battaglia et al] 
(2009j), where we omitted works directly addressed in the 
following. 
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1.3. REALM AF 

In this work, we present the code REALMAF (REAL MAgnetic 
Fields), which is aimed at measuring magnetic power spectra 
in galaxy clusters from maps of Faraday rotation measurements 
and apply this code to the Hydra A cluster. 

As will be detailed in the following sections, our method re- 
lies on these informations and assumptions: 

- the linearity of the measurement process (see Sect. 12.1b 

- the maps of rotation measures and their errors (see Sect. 14. Il l 

- knowledge of the geometry of the Faraday active volume, 
which also depends on the projection angle 6 of the jet and 
the cavity size (see Sect. 14.21 1 

- assuming that there is no polarised radio emission inside the 
Faraday active volume in front of the radio bubble (see Sect. 

- knowledge of the electron density profile Ug of the cluster 
atmosphere (see Sect. I4.2l i 

- assuming that the magnetic field stre ngth B follows the elec- 
tron density profile like B oc n" (e.g. lDolag etani200ll and 
see Sect.[3l1l 

- assuming Gaussian probability distributions for the magnetic 
fields and the measurement errors in the rotation measure 
map (see Sect. 12.2b 

- assuming statistical isotropy of B (see Sect. 12.3b 

- assuming solenoidal fields (see Sect. 12.3b 

- using a Bayesian approach to inference with a scale- 
independent prior for the power spectrum (see Sect. 13.1b . 

We use Hq = 70 km i ' Mpc ', ^ 0-3, Qa ^ 0.7 an d a red- 
shift of 0.0538 for Hydra A (Ide Vaucouleurs et al.l [T99l. which 
gives a linear to angular conversion factor 1" = 1.047kpc. To 
foUow the tradition in cosmic magnetism studies, we use cgs- 
units in this paper with the energy unit 1 erg - 10"^ Joule and 
magnetic field strength unit 1 G = 10 '^T. 

2. Magnetic field statistics in Faraday screens 

2.1. Measurement equation 

The rotation measure differs from the Faraday depth as defined 
by Eq.[T]only by the measurement noise n: 

RM(x^) = ,^(x^) + n(x^) 

= ^ I dz ap ne(xj_,z) Sj, Biix^,z) + «(xj_). (2) 

«xx (s.i) 

We assume the signal-to-noise ratio of the polarised flux to 
be sufliciently high for an accurate RM determination within 
an area A in the sky. The n-tuple of the RM values there then 
forms our data vector d = (RM(xj^))x^e4, from which all deduc- 
tions will be made. The RM data can be regarded as the result 
of a functional linear response operator R acting on the three- 
dimensional magnetic field configuration and producing a two- 
dimensional sky image of it, 

d^RB + n, (3) 

where the elements of the response matrix R are defined in 
Eq. |2] and the application of 7? on Z? should be read as the 
index-sum and line-of-sight integration given in the same equa- 
tion. The extended Faraday eff'ect in the intra-cluster medium in 
front of the Hydra A r adio galaxy was accurately measured by 
lTavlor&Perlevl(ll993h . The intra-cluster medium has a roughly 



spherical symmetric electron density ne{x.^,z) centred on the 
Hydra A cD galaxy. The radio jets of Hydra A inflate bubbles 
with relativistic, radio emitting plasma, which are visible as cav- 
ities in the X-ray emissi on of the surrounding hot cluster atmo- 
sphere (I Wise et al] 12007 ). On its way to the observer the po- 
larised radio emission of these bubbles experiences a rotation of 
its polarisation plane because of the Faraday effect of the mag- 
netised intra-cluster medium. 

2.2. Gaussian random fields 

Owing to the information lost in the projection, Eqs. |2] and [3] 
are not directly invertible. It is therefore impossible to calculate 
the line-of-sight magnetic field B^ from (f> and the known elec- 
tron density profile. However, we will see that assuming statis- 
tical isotropy of the magnetic fields and using a statistical ap- 
proach we are able to retrieve the amplitude spectrum of the 
magnetic field in Fourier space, the magnetic power spectrum 
cj(k) = (B(k) ■ B(k)), and the typical field strength B in position 
space and its autocorrelation length Ab, which are determined by 
a>(k). The expectation brackets should be understood as an aver- 
age over the ensemble of possible magnetic field configurations 
produced by typical intra-cluster dynamics. We will assume that 
the magnetic field statistics generated by those is Gaussian. 

Gaussian magnetic fields statistics is _nqt what numerical 
simulations of MHD turbulence find, see iWaelkens et alJ (l2009l) . 
However, they are the starting point of any analysis of correla- 
tion functions. If all we know statistically about the fields is their 
two-point correlation, the assumption which only expresses this 
knowledge is a Gaussian of which this correlation matrix is the 
covariance matrix. Any other distribution function would con- 
tain more information in a Shannon-Boltzmann sense. 

Before we can gain additional information on higher order 
statistics, we need to understand more the two point correlation 
function, and this is possible even with an approximate assump- 
tion on the statistics. Non-Gaussianities could modify the result 
of a power spectrum measurement based on the assumption of 
Gaussian fields. However, this modification can be expected to 
be of a modest nature, much smaller than the systematic uncer- 
tainties we deal with because of the imperfect modeling of the 
galaxy cluster's atmosphere structure and the geometry of the 
radio source that was used to probe for Faraday rotation. Our 
method is relatively robust to non-Gaussian modifications, be- 
cause the observational quantities it exploits, the product of RM 
values at different locations, are a linear transformation of the 
magnetic two-point correlation, irrespective of whether the un- 
derlying magnetic field statistics are Gaussian or not. 

Finally, many other methods also use Gaussia n random 
fields or assume thern im phcitly, e.g . iMurgia eta l. (200j); 
Bonafede et aP (l2010allbh : IvaccaetalJ (l2010h : iGovonietd] 
(2010). 

Using again the functional vector notation, we can write the 
Gaussian probability to find a magnetic field configuration B 
given a magnetic covariance matrix M - {B B^}, with elements 
M,/x,y) = <B,(x)B/y)): 



P(B\M) = g(B,M) = |27rMr5 exp|--B'^M"' By (4) 

with \M\ being the determinant of the magnetic covariance ma- 
trix. 
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2.3. Magnetic correlation tensor 

If there is a statistical translation invariance, M,y(x,y) depends 
only on the distance x - y because the absolute position can- 
not matter and only relative positions determine the two-point 
statistics. As a consequence, the magnetic covariance matrix be- 
comes diagonal in its Fourier space representation with respect 
to the k- vector: 



<B,(k) Bj(k')} = (Inf 6{k - k') M,/k). 



(5) 



If we deal with statistical isotropy, one moreover finds for 
solenoidal fields (V ■ B = 0) 

(jL)(k) I kikj\ ki 
M,/k) = ^ 1 - ^ + is^ji H(k), (6) 



(see e.g. lSubramanianl(ll999h . [EnBlin & Vod (l2003l) ) and there- 
fore 

aj(k) = Yj ^"(k)- (7) 

H{k) specifies the magnetic helicity, which does not play any 
role in this work, because it does not lead to any detectable sig- 
nature in Faraday rotation maps alone, as we will see shortly. 

2.4. Faraday statistics 

Assuming Gaussian statistics for turbulent magnetic fields has 
the pleasant property that the correlation function of the Faraday 
depth <p is also Gaussian: 



DBd(<p-RB) P(B) 



r DB r 

J llnMli J ' 



Dj] exp 1 2n i Tf{(p -RB)- - B'^M'^ B 



Dt] exp 1 27T i Tf(p - tj^R'^M R rj 



J^D?/ exp |2 



= \2nR'MR\-^ ex^\--^'{R'MRy 
= g{4>,R^MR). 



(8) 



Here, we used the Fourier representation of the Dirac delta func- 
tion (in functional spaces) as well as the well known prop- 
erty of Gaussian integrals J Ot]/ exp(-i/r^A"' + p-i]/) — 
\2nA\^l^ exp(y^A il2), which also holds for functional or path 
integrals denoted by The 0-covariance - {(f> (p^) - 

R^MR reads for our specific response and magnetic covariance 
matrices 

C^{xj_,x'J - al I dz I dz ne(x)ne(x') {B.(x)B^(x')} , 

M-(x,x') 

(9) 

where zi and Z2 denote the depth to the surface of the radio 
source at sky position and x^, respectively. If the spatial elec- 
tron distribution is constant, the magnetic fields are statistically 
homogeneous and isotropic and the observation samples a virtu- 
ally infinitely extended Faraday active volume, this can be sim- 
plified even more, and reads in Fourier space 



(^(k)0(k')> = (2jT)'6(k-k')C^(k) 
1 

2' 



C4k) - ^ alL.nl a>{k). 



(10) 
(11) 



where L, is the depth of the Faraday screen. Thus, clearly the 
magnetic energy, but not the helicity spectrum can be retrieved 
from Faraday rotation maps. In the more realistic case the geom- 
etry is less simple (finite observed volume, structured electron 
density, and magnetic field strength profiles, etc.). This formula 
has to be repl aced by one wh i ch tak es these effects into account, 
as derived in EnBlin & Vogd (I2003h . Because then the different 
Fourier modes couple, it can become simpler to again work in 
the position space, as we will do below (see Section [3T4l i. 



2.5. REALMAF and otiner power spectra estimators 

We model the magnetic power spectrum as a combination of 
spectral basis functions 



c^(k) = ^ SiWi(k), 



(12) 



where the spectral amplitudes form our spectral signal s = (s,). 
The number of bins, and therefore the number of free parame- 
ters of our inference problem, will be set by a trade-off between 
the aim to have high spectral resolution and to have small and 
uncorrected vertical error bars in the spectrum. 

A probability function for s given the data d, the poste- 
rior P(s\d), will be set up later in this work in Section 13.21 
Maximising the posterior with respect to the parameters in s 
provides an estimator for a>(k) expressed as amplitudes of rel- 
atively independent frequency bands. A simila r approach has 
alread y been applied to the Hydra A cluster by IVogt & EnBlinI 
(12005 '). Expressing the spectral basis functions as analogous ba- 
sis functions for the magnetic autocorrelation allows us to cal- 
culate the RM covariance matrix completely in position space, 
w hich makes two rough simplifications in the processing used 
bv lVogt&EnBiinl (l2005h unnecessary (see Section |5] for further 
discussion). Moreover we include noise effects into the analy- 
sis and model the Faraday screen more accurately using recently 
improved measure ments of the geom etry of the cavities in the 
Hydra A cluster bv lWise et all dioOTl) . 

Th e Hydra A cluster has been investigated by iLaing et aP 
( l2008h as well. They compared simulated rotation measure 
structure functions with the observed ones. A comparison with 
their results will be discussed in Section|5] 

Another approa ch for retriev ing magnetic power spectra is 
the one proposed bv lMurgia et alJ (.2004). Assuming a power law 
power spectrum, they simulated Faraday rotation maps by using 
the FARADAY code and compared the observed and synthetic pro- 
files of the dispersion of the RM as fit criterion. In order to find 
a good fit, iMurgia et alj (12004') assumed a power law to restrict 
the nu mber of degrees of fre edo m. This method was a lso applied 
e.g. bv lGovonietaP (l2006l) and lGuidetti et alJ (l2008l) . 

The REALMAF code presented in this work should alleviate 
several of the restrictions and assumptions of former approaches. 
Although REALMAF also requires the spectral range to be set up 
by hand, the posterior statistic calculated by it guides the choice. 
The power spectrum within the chosen range is relatively free 
to adopt any slope required by the data and permitted by the 
number of assumed spectral bands. 



' The algorithm described in Sect. in l3.3l mav not converge and may 
fail in cases where the variance of the spectrum is too large, or, equiva- 
lently, the inverse problem is too degenerate. For the our data this limits 
the spectral resolution to a maximum of about 7 bins. 



4 



Petr Kuchar and Torsten EnBlin: Magnetic power spectra from Faraday rotation maps 



3. Inference method 

3.1. Statistical inference and frequentist estimates 

In this work concepts from statistical inference will be applied 
to Faraday rotation data to measure magnetic field properties. 
Because this is traditionally done via frequentist estimates, we 
here briefly recall the conceptual differences in the methods. 

In a traditional frequentist data analysis the question is asked 
how likely the data would have emerged under certain model 
parameters, in the following called the signal. The parameters 
or signal values with high likelihood are then favoured over the 
ones with low likelihood. The likelihood is often estimated by 
comparing a statistical summary of the real data with that of sim- 
ulated mock data sets and by characterising their similarity and 
difference. In our case the Faraday rotation dispersion and em- 
pirical correlation functions are typically used as statistical sum- 
maries of the RM data. A comparison of data and model seems 
to be transparent, because it happens completely in data space. 
One can easily inspect the mock data by eye and compare them 
to the real one. This is one reason why frequentist methods are 
popular. Another one is their lower computational complexity, 
because they usually only require calculating the measurement 
process in a forward fashion. 

In statistical inference, a different question is addressed: 
How probable is a signal or a model given the data and other 
information sources? This is the reverse from the frequentist 
question and is actually the more relevant question scientifically. 
Answering it requires some knowledge or assumption on how 
probable the model was before taking the data. This prior proba- 
bility, in combination with the likelihood probability, permits us 
- using the Bayes theorem - to construct the posterior probabil- 
ity of a model given the data and any earlier knowledge, which 
then answers the question. Because this deduction takes the re- 
verse causal way of the measurement process, the computational 
complexity is usually higher than that of a comparable frequen- 
tist method. 

Thus, the statistical inference (or the Bayesian) approach is 
more abstract in trying to do deductions in the model or signal 
space. Its answer is a probability function in model space and 
not a single estimate. The mean and variance of the posterior 
are often used to characterise the solution. And no mock data 
and statistical summaries are usually generated, because the in- 
ference inserts the observed data into a likelihood function from 
which a backward deduction of the model is done. 

Frequentist and Bayesian results often agree. However, there 
are at least three reasons why they can differ, and we believe that 
the latter two matter here: 

1. Bayesian methods specify the prior assumptions explicitely, 
while those are often implicit in frequentist estimates. 
Because the prior has some impact on the result, different 
choices of priors will lead to different results. 

2. Bayesian methods are better suited to exploit the full infor- 
mation content of the data. The summary statistics used in 
frequentist analysis often looses a significant amount of the 
information imprinted in the data. In our case, we are inves- 
tigating the correlation between all pairs of pixels, but a typ- 
ical frequentist statistics might only look at the variance of 
the individual pixels or the correlation function from a sub- 
set of similar pixels, and therefore exploits much less of the 
data complexity. 

3. There are situations where no suitable summary statistics 
exist to construct an empirical likelihood. In our case, cal- 
culating an empirical RM correlation function to compare 



data with mock data only makes sense if the RM field can 
be assumed to result from a statistically homogeneous pro- 
cess. However, we are facing a very inhomogeneous geome- 
try here because of the highly structured cluster atmosphere, 
with very few comparable lines of sight. How to obtain a 
meaningful two-point correlation function from the data in 
this situation is not clear to us. 

Because of the last point, the Bayesian inference performed 
in this work does not provide common quantities of frequen- 
tist analysis like mock data, observed RM correlation function, 
and the like. The RM data are confronted with the theoretical 
covariance matrix that is expected for a given power spectrum 
of the magnetic field after a model of the inhomogeneous clus- 
ter atmosphere is factored out. The power spectrum would, in a 
statistically homogeneous setting, describe an easily observable 
RM correlation function, but not in our inhomogeneous case. 

The Bayesian approach also extracts statistical information 
from poorly sampled quantities, while providing the correct 
large uncertainty resulting from the undersampling. For exam- 
ple in our case the magnetic modes on scales of or larger than the 
observed RM map are clearly undersampled. The largest modes 
appear as a single number in the data, the mean RM value of 
the map. Because it is a single number, it is hard to invert this 
into precise spectral information on large-scale magnetic fluc- 
tuations. Nevertheless, there is some useful information in this 
number, especially if it is far from zero. In this case, the total 
magnetic spectral power on large spatial scales should be on 
a level which makes the observed mean RM value sufficiently 
likely. The use of a Jeffreys prior on all spectral amplitudes en- 
forces in our application that the inversion of poorly constrained 
modes is made in a very conservative fashion. With frequentist 
methods the usage of the information in such poorly sampled 
modes is more difficult, because no empirical statistics can be 
built up from a single observation. 

3.2. Statistical basics 

Our aim is to retrieve from the rotation measures (data d) the 
magnetic power spectrum represented by the spectral amplitudes 
(signal s). The Bayes theorem expresses the posterior P(s\d) (the 
conditional probability of the signal given the data) in terms of 
the likelihood P(d\s) in the following way: 



P(s\d) = 



P{d\s) Pjs) 
Pid) ■ 



(13) 



Because P{d) is a constant with respect to s, the most probable 
s given d can be retrieved by maximising the joint probability, 
which consists of the likelihood P{d\s) multiplied with the Prior 
P{s). These will be specified now for our problem. 

In our case the signal is the ensemble of spectral amplitudes 
s - {si)i. We adopt Jeffreys prior for each spectral amplitude, 
so that P{s) oc Yli sj^ . Jeffreys prior is a uniform distribution on 
a logarithmic scale, expressing that we are beforehand not even 
sure about the magnitude of the magnetic spectrum. 

Our likelihood Pid\s) has to incorporate two stochastic pro- 
cesses. The first is the realisation of the magnetic field con- 
figuration according to Eq. H) Because s fully determines the 
here relevant symmetric part of the magnetic correlation tensor, 
M - M(a),H) - M(Y,i SiCJi,H), the statistics of the Faraday 
depth is determined by it as well, and it is also Gaussian ac- 
cording to Eq. [8] If the noise is also Gaussian, then the data are 
Gaussian with a covariance matrix C, with entries Cij = {dt dj) 
for the correlation of RM data at locations x^., and Xj^ y. Using 
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Fig. 1: Left: the magnetic autocorrelation tij(r) in position space for a typical magnetic power spectrum reconstruction for the cool 
core in the Hydra A cluster using seven spectral bins. The negative part of w(r) is dotted. Right: the counterpart in Fourier space 
a)(k). The parameterisation can adopt any shape permitted by the seven basis functions seen as wiggles in the power spectrum. The 
parameters 6 = 45° and a - I were assumed, with 6 described in Section l4r2] and a introduced in Eq.fTTl 



vector and matrix notation for the data vector d - 
covariance matrix C, we can write the Gaussian as 

P{d\s) = g(d,Cis)) = \2nC\-^- g-J"'^'^"'''. (14) 

An expression for C as a function of our spectral parameter s is 
needed. In the view of the linearity of Eq. |2l the observed data 
RM can be expressed in terms of the magnetic field B (vector 
over physical space), a response matrix R (projecting from phys- 
ical to data space) and noise n (in data space) as. d = R B + n. 
From this we find 

C(s) = {iRB + n)iRB + nf) = R{B B^) R^ + {nn^) 

= RMR^ + N = C4s)+N. 
Thus to any derived Faraday rotation covariance matrix - 
{<p<p^) an estimate for the noise covariance matrix - {nn^) 
needs to be added. 



3.3. Maximising tine posterior 

The power spectrum inference within the maximum-a-posteriori 
approach can only be conducted numerically. Here, we max- 
imise the logarithm of the posterior, because there are con- 
venient analytical expressions for the gradient and Hessian. 
Multiplication of the likelihood with the prior to get the pos- 
terior is an addition in logarithmic space, which also simplifies 
matters. Indices /, j refer to the degrees of freedom of the param- 
eterisation of s. We find 



n and its 3.4. Tine covariance matrix 



d\n[P{s\d)] 

dsi 



iTr 

2 



K-c)(c-.gc-) 



(15) 



and if s contributes only linearly to C: 



d'^\n[P(s\d)] 

dsjdsi 



ds 



dsi 



(16) 



-Tr 



,.2 ■ 



The most probable power spectrum is then found iteratively. We 
use a second-order Newton method to solve for it, with 

d^\n[P{s\d)]V^ d\n[P{s\d)] 



dsjdsi j dsi 
where repeated indices are assumed to be summed over 



An expression for the covariance matrix - {(p(xj_) 0(x^)) of 
the rotation measure map needs to be derived in the following 
form: 

C^{x^,x'J = al < £ dzn,{x)B,(x) £ dz' n,(x')B,(x') ). 

The electron density distribution «j,(x) needs to be obtained from 
independent data, e.g. from X-ray observations. T he magnetic 
field e nergy is assumed to follow ne(x) according to lDolag et al.l 
(I2002h : 



<B^(x)) = B^(«,(x)/no) 



(17) 



The following approach is based on statistical homogeneity of 
the magnetic field. Therefore we re-scale (B^(x)) with a rescale- 
ing function h, h{x) - {B^{xyf -^ I Bq - (ndx) / no)" . This ensures 
{B\x)} - {B^{x))lh^{x) - Bq to be independent of position. If 
we also assume statistical isotropy, we obtain the scaling of the 
required z-component of B: {B\{x))IBI^ = {B^(x))/bI. Thus, 

/-*oc 

= al I dzne{xj_,z)h{xj_,z) 

dz riAx'^, z')h(x'^, z'){B,(x^, z)B,(x'^,z)) 
alnl J dzf(x^,z) J c/r, /(x^, r, + z) M,,(r), (18) 



r 



with r = (xx -x^, Tj), where we introduced the window function 
/(x), which is nonzero inside the Faraday screen: 



fix) = h(x) = 

no \ no 



(19) 



M~.{r) - {B,{x)B.{x+r)} is the zz-part of the (rescaled) magnetic 
autocorrelation tensor, which we try to measure. 

3.5. Modeiiing tine magnetic autocorreiation 

M.~{r) can be split up into a longitudinal part M]^(r) and a cor- 
responding perpendicular part Mj^{r), which depend only on the 



6 



Petr Kuchar and Torsten EnBlin: Magnetic power spectra from Faraday rotation maps 



absolu te value of r (see ISubramanianI ( 1 19991) or lEnBlin & Vogtl 
(l2003h V 

r? rl 
M,-Sr) = ML(r) ^ + M^{r) ^ (20) 

The rescaled magnetic autocorrelation function 
w(r) = {-S(x) ■ B(x + r)) is given by 

w(r) = 2 M,v(r) = lM^{r) + ML(r). (21) 

The information of the data can be combined into a one- 
dimensional correlation function and therefore does not allow us 
to simultaneously determine two independent functions. Thus a 
relation between Mn and Ml is necessary. There are two suitable 
possibilities how to connect them. 

a) Full isotropy assumption, which additionally assumes 
isotropy at any location in Fourier space. This yields 

MN(r) = ML(r). (22) 

b) Considering the solenoidal character of the magnetic field by 
V ■ B = 0, which yields 

MN(r) = l^(r2ML(r)). (23) 
Lr ar 

The V • B = condition leads to Eq.|6]and is more physical, but 
also more complicated to handle. Both assumptions have been 
implemented in REALMAF. However, the most computationally 
efficient way is to use the full isotropy assumption and apply the 
divergence relation in Fourier space, which results in a simple 
multiplication factor, see Section lTSl 

Connecting of Ml and Mn allows us to parameterise M,,(r) 
using 

M,.(r) = ^ SiM^ir), (24) 

where M^lV) is further split into wf^{r) and M*j'(r) described 
by Eq.|20] Because of the linearity of the Fourier transformation, 
the linear coefficients s -, in real space also appear as linear coef- 
ficients in the analogous representation in Fourier space. Thus, 
although we are calculating the covariance matrix completely in 
real space, the power spectrum can be obtained directly: 

oj(k) = Yj (2<'(^) + <''(^)) = 2 ^/ 

. The basis functions together with the retrieved spectral ampli- 
tudes Si can be translated into w(r) (see Fig. [T| a) and the cor- 
responding (L>{k) (see Fig. [T| b). The used model functions for 
M*j' and M® were chosen to transform into relatively well sepa- 
rated and well localised frequency bands Wi{k) in Fourier space. 
Appendix|A]shows the used model functions and describes how 
they were derived. 

3.6. The divergence relation in Fourier space 

The zz-part of the magnetic autocorrelation tensor can be ex- 
pressed in Fourier space analogously to Eq. |20] 

^ k^ 

M,,(k) = Mdk)^ + Mn{k)-^. (25) 

Applying V ■ B = gives then 

M^(fe) = ^ a){k) = 2Mn{k), (26) 



see lEnBlin & Vogtl (l2003l) . The M^ik) from Eq.|25]is in general 
not the Fourier transformation of ML(r) used in Eq.|20] To create 
a connection to the real-space parameterisation M,,(r), we mul- 
tiply aiik) by 3/2, which is analogous to an addition of Miik) 
with Mi(k) - Mf^ik). This makes My,(Vi) a spherically symmet- 
ric function, which forces that its Fourier back-transform M^{r) 
must be spherically symmetric as well. Then, 

ML(r) = MN(r) ^ M^k) = Mn{k\ (27) 

which is the full isotropy assumption defined in Eq. l22l 
Therefore, when assuming full isotropy, the divergence condi- 
tion can be imposed on co{k) to a good approximation via multi- 
plying it with the factor 2/3. This was verified as well by numeri- 
cal experimentsQ Using the full isotropy assumption reduces the 
processing time by up to a factor 4 and allows a simpler spher- 
ically symmetric model for M,,. In practice the power spectrum 
is multiplied with a factor 2/3 and the resulting magnetic field 
strength by V2/3. 

Finally, the V ■ B - condition, which is used by 
,Vogt&Enmin (2005l) as well, is not exact because of 

V B^hV B + B Vh^B Vh, 

with the rescaleing function h defined in section [3741 If B varies 
on scales much smaller than the ones on which h varies, this as- 
sumption is sufficiently close to the condition V ■ B = 0. Another 
indication that the resulting errors are negligibl small is that a 
totally different assumption on the statistical distribution of the 
magnetic fields, like full isotropy, which violates V ■ B = by 
definition, produces only a shift by a factor of 1.5, as shown 
above. 

4. Application to the l-lydra A data 

4.1. Tlie data 

High-resolution polarisation measurements of the radio lobes of 
Hydra A at different wavelengths, wi th the VLA suitable fo r 
RM map making, were performed by iTavlor &Perlevl(ll993h . 
High-fidelity RM maps were constructed from the se using the 
PACERMAN code (Dolag et al. 20051: PVbgt et al.'2005). PACERMAN 
uses error weighting and also non-local information to solve the 
nTT-ambiguity and therefore produces more reliable Faraday ro- 
tation maps than algorithms that do standard fi tting of Eq . |2lpixe l 
by pixel. The resulting RM map was tested bv lVogt et al.l d2005l) 
for the presence of map-making artifacts. These artifacts can be 
statistically detected using the anti-correlation of inferred RM 
and position angle values imprinted by misfits, as first shown by 
lEnBlin et al] (l2003b . The tests show that the current RM map of 
the south lobe of Hydra A is prone to n;T-ambiguities and other 
artifacts, even if generated with PACERMAN. However, the north 
lobe, which contains lower RM values an d is therefore pointing 
towards us, is sufficiently clean. Therefore lVogt & EnBlinl(l2005h 
analysed only the north lobe of Hydra A, and our application of 
REALMAF is also restricted to this area. 

The resolution of the RM map is too high to take every single 
pixel into account in our analysis. A computational bottleneck is 
the double integral in the calculation of the model covariance 
matrix followed by matrix multiplications described in Section 
13.31 To reduce the amount of required computational resources, 
we merge areas in the map to bigger pixels by combining the 

- Actually this relation was discovered in numerical experiments 
first. 
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Fig. 3: Left: Radial profile of the original rotation measure data and the merged values shown in Fig. |2l To facilitate reading of 
the plot, the 1000 merged data points are combined within (variably sized) radial bins. Right: Comparison of the merged profile 
with the theoretical radial profile of the ensemble average resulting from the analysis. The retrieved profile was obtained from the 
trace of the covariance matrix described in Eq. [T8]and implicitly emerges from the modelled magnetic power spectrum shown in 
Fig. [8] and the assumed geometrical profiles of magnetic energy and electron density. The former determines the amplitude of the 
theoretical profile, whereas the latter determines the detailed radial decline. This theoretical curve was not obtained by a fit of the 
radial profile, which is shown here for illustration only. Instead, our inference is based on examining the cross-correlation of all 
pairs of points, which contain much more information compared with the autocorrelation of individual points displayed here. We 
note that the variations seen in the data with respect to the theoretical mean curve are a natural consequence of the data being only 
a single realisation of a random process, whereas the theoretical curve represents the ensemble average of a large number of these 
reaUsations. 




Fig. 2: Merged data map of Hydra A with 1000 pixels remaining. 
The RM values range between -1700 rad/m^ (dark) and 2800 
rad/m^ (light). The map has a size of about 40 kpc. Each square 
represents one value used in the analysis. 



rotation measure values, the errors, and the pixel position in an 
inverse noise-weighted fashion. The idea is to average the data 
in quadratic areas chosen so that they have roughly the same 
noise level. Therefore the outer regions with lower radio surface 
brightness and consequently larger RM errors are more heavily 
averaged. For details on the algorithm combining the pixels, the 
reader is referred to Vogt & EnBlin (2005). 

These merged data points then form the data vector d used 
in section [3T2I Figure |2] visualises the merged map by showing 



the averaged values in colour within the square areas from which 
they were averaged. The area outside the squares does not con- 
tain any used RM values. In Fig.[3]the radial profile of of RM^ 
is shown and is compared with the ensemble averag ed (RM^) 
predicted by the deduced magnetic spectrum and assumed radial 
profile of the electron and magnetic field energy densities. From 
this figure it can be seen that the largest separation between our 
pixels is r^ax = 33 kpc, whereas the smallest r^i^ = 0.073 kpc. 
A first guess of the optimal spectral range can be obtained based 
on these maximal and minimal sizes of visible structures calcu- 
lated using fcmin - Tr/^max and A:max = ^/'"min The final result in 
Fig.Quses roughly these spectral boundaries. 

Spatial averaging of RM data is not properly modelled in 
Eq. [18] In principle, it could be included into this formula and 
the generic inference formalism, owing to the linearity of aver- 
aging. However, we have to perform the double integral in Eq. 
[TS] numerically for any pair of data pixels. Adding four more 
integration dimensions would exceed our computational possi- 
bilities. 

We accept the spectral error resulting from the data averag- 
ing operation because pixels that resulted from averaging over 
large regions are automatically only contributing information 
about scales above their size. Any baseline formed by them with 
other pixels will be of the order of the pixel size or larger. The 
product of the RM values of a pair of pixels is compared with 
the expected data correlation over this baseline distance. Thus 
these pixels determine the inferred spectrum mostly on scales on 
which their internal RM substructure, which was averaged out, 
does not matter. 

The silhouette of the radio lobe used to probe the Faraday 
screen breaks the translational invariance of our problem, in ad- 
dition to the effect of the electron density profile. However, this 
does not impact directly our data covariance matrix C, because 
this only contains the expected (averaged over an ensemble of 
data realisations) correlation of individual pairs of positions in 
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the RM map. The silhouette restricts the available positions and 
therefore the number of baselines, which means that on larger 
scales, we have less information. However, the Bayesian for- 
malism is informed about this and just increases the uncertainty 
margins for poorly constrained data modes. 

We also require to specify the noise covariance matrix. 
Assuming the noise at different pixels to be uncorrelated, the 
noise matrix is diagonal. The diagonal elements can be approxi- 
mated by the square of the RM-deviation of each pixel given by 
the error map produced by PACERMAN. 

The used rotation measure map has a mean value of about 
600rad/m^. The galactic rotation measure at the position of 
Hy dra A is about 1 percent of this and therefore negligible, 
see iJohnston-HoUitt et al.l (l2004t) . Therefore we assume that the 
mean value is caused by cosmic variance which is in turn pro- 
cessing a too small sample of the cluster. This mean value is 
correctly interpreted by our method in that it adds power to the 
small k part of the magnetic power spectrum, compared with the 
case in which we subtract the mean value from the data d. 
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Fig. 4: Histogram of the minimal distance 7?„„„ from the cluster 
centre for 800 data points depending on the assumed projection 
angle 9. Magnetic fields at r < /?„„■„ are invisible to the method. 



4.2. Window function 

The window function / defined in Eq. [19] requires the elec- 
tron density distribution ndx). Hydra A is a cooling-core cluster 
Therefore we assume a spherically symmetric double beta pro- 
file defined as 
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(28) 



with «ei = 0.056 cm ^, n^^ - 0.0063 c m ^, rd = 33.3 kpc , 
r,2 = 169 kpc ar id /? ^ 0.766 as fi tted by Vogt & E nBlinI (l2005h 
using data from iMohret all (119991) . Various works propose that 
the observed radio lobes of active galactic nucl ei, like Hydra A. 



form c a vities in t he clu st er atmosphere, se e e.g.lBohringer et al 



T^, lEnBhnl (fl99l iMcNamara et al.1 (12000^ 



Fabian et al 



2000"), and 'David et alj (2001). In these cavities the polarised 
radio light is produced, but «(.(x) is presumably so low that no 
significant Faraday rotation is expected to occur there. The fact 
nj^x) - for any x € cavity can be taken into account by adjust- 
ing the integration boundaries in Eq.[T8] In order to calculate the 
important lower integration boundary, where the Faraday active 
area starts, we need to know the geometry of the cavity. The cav- 
ity in our used map-section can be approximated by an ellipse 
with the projected distance 24.9 kpc from cluster centre, pro- 
jected s emimajor axis 20. 5 kpc, and semiminor axis with 12.4 
kpc, as Wise et all (l2007l) propose. Additionally the Hydra A 
north lobe is tilted by an angle Q towards the observer. Depending 
on the innermost radius of the probed Faraday active area 
shifts, as can be seen in Fig. |4] There it becomes apparent that 
the sensitivity of our analysis starts at radii of about 20 to 30 kpc 
from the cluster centre depending on Q. 8 is not well determine d 
and its estimates range from about 30° (lTavlor&Perlevl(ll993h 
if the c avity eff'ect is add ed) to more than 45° dLaing et al.ll2008l) 
to 60" (IWise et al.ll2007h . 

4.3. Measuring fhe magnetic profile 

The index a of the magnetic field can be retrieved from the data. 
In contrast to the power spectrunfl which is built dynamically 
by our algorithm, a has to be set up for every processing run. 




and the directly obtained central magnetic field strength Bo 



0.35 0.4 0.45 0.5 0.55 0.6 0.65 
kmin [kPC"^] 

Fig. 5: Sampling of k„,i„ for given a = 1 



Maximising the posterior with respect to a provides indications 
on most probable configurations. One drawback is that a also 
depends on the low-frequency cut-off k,„i„ (position of the top of 
the first spectral bin) chosen; meaning a two-dimensional prob- 
lem needs to be optimised. The result we find is a = 1 with 
kmin = 0.42 kpc"'. Because this result only depends on the large 
scales, a rather low-resolution study with 300 data points was 
sufficient and permitted a fast investigation of the 2-d parameter 
space. Sampling of k,„i„ for three diff'erent projection angles 6 
can be found in Fig. [5] Cut-offs with k < 0.42 kpc ' coiTespond 
to a > 1 and vice versa. Smaller 6 implies a smaller physical 
distance between the data points, which results in a small shift 
of kmin- However, it is negligible and inside the statistical vari- 
ance of about 0.05 kpc Fig. |6] shows a consistency check by 
sampling a for the given fc„,„, = 0.42 kpc"' . The curvature in Fig. 
|6] gives roughly a 1-sigma variance for a of about Aa = +0.3, 
while assuming Gaussianity around the peak. When we also con- 
sider the impact of the variance of kmin, we find a - I + 0.5. 

Surprisingly the posterior maximum is rising for rising 6 
with no apparent limit. We suppose that this has nothing to do 
with the probability of the physical model and therefore 6 can- 
not be restricted using the data of the Hydra A north lobe only. 
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Fig. 6: Sampling of a for given A:,„,„ - 0.42 kpc 
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Fig. 8: Inferred power spectrum for three assumed projection an- 
gles and comparison with a Kolmogorov slope (dotted). Note 
that for k > 10 kpc"' beam effects reduce the amount of power 
and are probably responsible for the steepening. 
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Sso [mG] 


Ab [kpc] 


spectral index 


30" 


21 ± 1 


9.3 


5.0 ± 1.0 


1.70 ±0.14 


45° 


36 ±2 


16 


5.2 ± 1.1 


1.73 ±0.13 


60° 


85 ±5 


37 


5.3 ± 1.2 


1.85 ±0.14 



Table 1 : Magnetic field characteristics for most probable a - I 
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Fig. 7: Power spectrum with an optimal, a too small and a too 
large cut-off ^,„„, 



Indeed, is a parameter, that only affects the integration bounds 
of Eq.[T8] whereas the power spectrum, the magnetic scaling a, 
and the cut-offs affect the integrand. 

4.4. Power spectrum and magnetic characteristics 

The shape of the power spectrum depends on the spectral cut- 
offs k„i„ and k,„a_^. We retrieved the most probable k,„i„ to be 0.42 
kpc"' in Section l43] Obviously, the most probable configuration 
is a power law with no large-scale turnover visible within the 
probed spatial scales, as can be seen in Fig. |7] The decline for 
smaller k„,i„ may be caused by the limited rotation measure map 
dimensions and the down-pull of thereby unconstrained ampli- 
tudes by Jeffreys prior A larger A:„„„ will remove necessary de- 
grees of freedom of the model at the large scales, which forces 
the algorithm to deposit additional power at the neighbouring 
scale. The high-frequency cut-off is limited by the resolution of 
the map and instrumental beam. A beam with spatial scales of 
0.32 kpc (FWHM) gives k^ax = 10 kpc"' . Below, we increase the 
resolution by increasing the number of data points to 1000 and 
extend the spectrum to a maximal possible range. We made three 
processing runs with a = 1 and three possible projection angles, 
see Fig. |8] and Table [1] To investigate the dependence of the re- 



sults on a we made two more runs with a - 0.5 and a - 1 .5, see 
Table |2] Formulas for the calculation of the magnetic character- 
istics are described in Appendix |B] The given statistical uncer- 
tainties refer to static a and 6. These errors are propagated from 
the error bars of the power spectrum. 

The power law ranges from scales of 0.3 kpc (beam FWHM) 
to about 8 kpc. The fall-off in power at small scales is caused 
by instrumental beam effects. However, the larger scales are 
not greatly affected by the beam. To measure the spectral in- 
dex we ignore the first spectral bin, which is possibly affected 
by map boundary limits as well as the last two beam affected 
bins. The dotted line in Fig.|8]represents a Kolmogorov slope of 
e{k) oc A:"'-^^ compared with the plot using 6 = 45°. We achieve 
a perfect fit for the central part of the retrieved power spectrum, 
which is also inside the larger statistical variance of the parts at 
smaller k. There is a noticeable dependence of the spectral index 
on a with minimal and maximal slopes of 1.5 and 2.1, which 
should be regarded as improbable worst cases. 6 has a big im- 
pact on the resulting magnetic field strengths. Assuming a larger 
6 implies a shorter interaction path with the magnetic field and 
therefore greater necessary magnetic field strengths. The derived 
Bo in the cluster centre assumes that the magnetic profile given 
by Eq.[T7]can be extrapolated with a constant a into the cluster 
centre. However, the sensitivity starts at regions from about 30 
kpc from the cluster centre, see Fig. |4] and a flattening of the 
magnetic profile at smaller radii is well possible. Bgo gives the 
directly probed magnetic field strengths at 50 kpc distance. 

The impact of the instrumental beam on the calculated mag- 
netic field strength is low. A continuation of the power law to the 
small scales would give an additional 0.2 fiG at most. 
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a 


So UiG] 


Bso [pG] 




[kpc] 


spectral index 


0.5 


18 ± 1 


12 


4.8 


± 1.1 


1.56 ± 0.14 


1.5 


70 ±4 


21 


5.5 


± 1.2 


2.06 ± 0.14 



Table 2: Magnetic field characteristics for — 45° and maximal 
and minimal possible a 



author 


a 


Bq in i^G 


spectral index 


large-scale turnover 


Kuchar 


1.0 


36 


1.73 


no 


Laing 


0.25 


19 


0.8 


no 


Vogt 


0.5 


7.3 


1.67 


yes 



Table 3: Recent results of Hydra A for 6 — 45 



5. Comparison with other methods 

The magneti c power spect r um o f Hydra A was previousl y 
measured by Vogt & EniBliiil (l2005h and by iLaing et all (l2008h . 
Both find significantly lower central magnetic field strengths 
than we do. This is because of a much flatter magnetic field 
profile in ferred characterised by a smaller index a. Moreover 
IVogt & E nfilin ( 2005) found a bump, a cut-off at spatial scales of 
about 2 kpc and more, but the spectral index for small scales is 
still similar to our current results. However, Laing finds a much 
flatter spectral index of about 0.8. Table |3] compares the results 

for a projection angle of 6* = 45". 

To partly reconstruct the results of lVogt & EnBlinI (l2005h . we 
applied the two simplifications they have used in the calculation 
of the covariance matrix in Eq.fTSl 

1 . extending the second lower integration bound to -oo 

2. merging the two window functions / assuming the mag- 
netic autocorrelation changes on much smaller scales than 
the window function. 

Then we compared the results with the Fourier method of 
Vogt, but adding noise effects, the cavity structure, and a 
more advanced triangular parameterisation of the power spec- 
trum, see Fig. |9] We reconstructed the spectral bump found by 
IVogt & E nfilin (2005) and found a perfect match between both 
methods. Figure [10] shows the difference to our method without 
these simplifications. 

Consistent with IVogt & EnfihnI (|2005), the most likely a is 
shifted to much lower values if one adopts these simplifications, 
as can be seen in Fig.[TT] However, instead of a falling plateau of 
the posterior between a = 0.1 and a - 0.8 and a cut-off for the 
unphysical range c < 0.1 we again find a Gaussian distribution, 
but with a maximu m at 0.1. This very low a value indicates that 
the bump found by I Vogt & EnfilinI (l2005h may be explained by 
a power leak on large scales produced by an implicit increase of 
the large-scale contribution of the window fu nction. 

One of the major simplifications used bv lLaing et al.l (l2008l) 
is a linear and position-independent relation between the RM 
power spectrum and the magnetic autocorrelatiorQ, which is 
only valid for a deep and statistically homogeneous Faraday 
screen. Both assumptions are not very well met for the Hydra 
A cluster If we simulate a constant window function by setting 
a = - 1 (for which the effective window function, as given by 
Eq. [19] becomes flat) and use symmetric integration bounds in 
Eq. [18] we get a spectral index of about 1.1. For a - 0.25 we 
get a central magnetic field strength of 12.9 juG, which is much 
closer to their results. 

* applied bv lVogt & Enijlid (l2003h as well 




1e-16 



1e-17 I — — ■ — ' — ' — —I 

0.1 1 10 

k[kpc"''] 

Fig^9: Reconstruction of the results of the Fourier space method 
bylVbgt & Enfilin (2005) by adding the removed simplifications 
for a - 0.5 and 6 = 45". This plot illustrates the continuous 
result taking the model for the power spectrum into account. 
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Fig. 10: Power spectru m using the REALMAF and the method by 
Vog t & EnfihnI (l2005h for a = 0.5 and 6 = 45°. According to 
Fig.|9]the Fourier method is equivalent to explicitly adding sim- 
phfications to the real space method. 



The remaining difference to the spectral index of 0.8 and 
field strength of 19 juG m ay be partly caused by the beam 
correction applied by L aing et alj ( l2008i) . which lifts the high- 
frequency part of the spectrum, which is possibly contaminated 
with noise residuals from the RM-map-making algorithms. An 
implementation of a beam correction as a convolution in real 
space is not feasible in REALMAF, because it would be computa- 
tionally too expensive and an implementation in Fourier space is 
in general only possible with simplifying assumptions. However, 
the exact frequency-dependent instrumental beam is unknown 
anyway, and we believe its effect can be neglected at scales larger 
than the beam size. 



6. Conclusions 

6.1. Summary 

We developed a new method to retrieve magnetic power spec- 
tra from Faraday rotation maps implemented in the REALMAF 
(REAL MAgnetic Fields) code, where we modelled the mag- 
netic autocorrelation function in real space. 
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small scales. We published the REALMAF code for general scien- 
tific usage under an open source licens^ 

Upcoming telescopes such as EVLA, LOFAR, LWA, 
ASKAP, and SKA will certainly provide the necessary high- 
fidelity Faraday rotation data of many more galaxy clusters for 
further studies with REALMAF. This may help to figure out the 
origin of cluster magnetic fields and the role they play in the 
metabolism of the intra-cluster medium. 
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in IVoet & EnBIinI j2005l) available for insight, thanks to Daria Guidetti and 
Robert Laing for discussions and to an anonymous referee for constructive com- 
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Forschergruppe 1254 "Magnetisation of Interstellar and Intergalactic Media: The 
Prospects of Low-Frequency Radio Observations". 



Fig. 11: S ampling of a with and without simplifications used in 
IVogt & EnfiUn (.2005.) for 6 = 45" and k,„i„ = 0.3 kpc"' 



Furthermore we introduced and verified a way to take the 
solenoidal character of the magnetic fields into account by using 
a simple full spherical model for the magnetic autocorrelation 
and multip lying with a proportio nality factor In contrast to the 
method bv lVogt & EnBiinl(l2005l) with a model in Fourier space, 
this permits us to alleviate some previously required simplifi- 
cations in the processing. However, we are able to reconstruct 
their results by synthetically adding their simplifications to our 
method. 

When we apply REALMAF to the north lobe of the Hydra 
A galaxy cluster, we find a power law power spectrum in the 
spatial range of 0.3 kpc to 8 kpc, with no visibl e large-scale 
turnover within this range, as reported by I Vogt AEn61in(2005|). 
The magnetic density profile, which we modelled with an index 
a applied to the electron density profile, seems to follow the elec- 
tron density profile with ff = 1 ± 0.5. The somewhat large uncer- 
tainty on a also aff'ects the spectral index of the power law, which 
is then in the range of 1.5 to 2.1. The exact projection angle of 
the used radio lobe Hydra A north is not well restricted, but it 
only seriously aff'ects the absolute value of the magnetic energy 
density in the cluster and not the spectral slope. Depending on 
the angle we obtain field strengths ranging from 20 fiG up to 80 
jU G in the cluster centre. However, these are extrapolations, be- 
cause the measurement is only sensitive at a distance of about 
50 kpc from cluster centre, where we find fields between 10 //G 
and 30 //G. The extrapolation of the model into the centre may 
be questionable. 

Theories descri bing the generation of magnetic fields 
by turbulence like Subramanian et al. (2006), En Blin & Vogtl 
(2006) and Schekochihin & Cowley (2006) forecast a large- 
scale turnover, which might be visible inside our spectral range. 
Compared with IVogt & EnBlinI (l2005h we found no large-scale 
turnover within the spectral range up to 8 kpc. A turnover at 
larger scales is, however, not excluded by our analysis. 

6.2. Outlook 

The REALMAF code will permit us to analyse the Faraday rota- 
tion data of other galaxy clusters. This will hopefully allow us to 
study magnetic field spectra over a large range of spatial scales 
and in a variety of cluster environments. A usable map of Hydra 
A south would further increase the probed range at large scales 
and a map with a higher resolution would increase the range at 



Appendix A: Model functions 

Here we describe the construction of our spectral model func- 
tions. We start with an ansatz for M]^(r) for the case V ■ B = 
and use a Gaussian multiplied with a cosine. The corresponding 
MN(r) is obtained by Equation|23] Thus we adopt 

1 1 

Mj\r) -——b^ exp[ -] cos[/7,r] , and 

V2^ 2^2 

M«(r) ^ f'^ exp[-;^]((2fl2 _ r2)cos [bA ^^'^^ 

+a^ibirs,m [bir]). 

These functions are Fourier transformed analytically using 

r°° sin F/trl 

M{k) = An r^M(r) —dr, 

Jo kr 

see lEnBlin& Vogl(l2003l) . We bind the free parameter a, as 




because in this case (J-' \k) = 2M'^\k) + Ml\k) is fully positive 
and concentrated around a maximum at A; = 3 bi. This allows us 
to put the power spectrum together by relatively well separated 
and located frequency bands, as Fig. IA.2I illustrates. Fig. lA.ll 
shows the model functions in position space. Especially Mn(?") 
has a strong negative part, because any closed field line has to 
come somewhere back through the plane perpendicular to it. 

If we assume full isotropy, M.~(r) becomes spherically sym- 
metric. In this case a suitable model function can be obtained by 
Eq.IAH 

M»(r) = io.®(r) = i(M«(r) + 2M«(r)). (A.2) 
This implicitly assumes MN(r) - Mi^{r). Thus we get 

a) using V ■ Z? = 

in position space: 

b^ 8 
Ml\r) = —^exp[--b^r^]cos[bir] 
yTjT " 

M®(r) = — ^ exp[-lbjr\2(9 - Sbjr^)cos[b,r] 
18 yljT " 

-9fe,r sin[fe,r]), 

(A.3) 

' ,http : // sourcef orge . net/pro j ects/r ealmafT] 
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Fig. A.l: model functions mJ^'C") and M^^(r) with a,- - | and 




Fig. A.2: Model functions M^^'(k), M^'ik) and Lo^'\k) in Fourier 
space with a, = | and bi - 1 



in Fourier space: 



(,) 27 9(fc,+^)% 

1024b^k ^ 32b] 

[l6b] + 9b] k + 18bik^ + 9k^ + 9P 

9k \ 
+ exp[— ](-16fef + 9b]k-\%bik^ + 9k^)\ 



(A.6) 



Appendix B: Calculation of magnetic the field 
characteristics 

Here we provide the used formulae for the magnetic field char- 
acteristics. 

a) Magnetic field strengtli 

The magnetic energy density can be calculated by integrating 
over the power spectrum or analogously taking the autocorre- 
lation function in real space at the origin w(0). The magnetic 
field strength in the centre of the cluster, where f - I and 
therefore B - B'k then 



b) Magnetic autocorrelation lengtli 

Xoo 
-co ^'"'^('") 



w(0) 



(B.l) 



(B.2) 



c) 1 -dimensional power spectrum 

Transformation from the three-dimensional power spec- 
trum cojk) to the one-dimensional e{k) is done using 
(lEnBhn & Vogtll2Q03h 



e(k) 



k^u){k) 
2{2nf ' 



(B.3) 



d) Spectral index using a Bayesian approach 

The spectral index can be calculated from the slope of the 
power spectrum in logarithmic scale. The fitting formula is: 



^ ' 64 



27 



(k-bi ^ 9(bi-k)\ 



k + bi ^ 9(bi + k)\^ 
+ —, — exp[ —5 — ] 



27 



32^7/ 

(, 9(b^-kf m-kf 

(—TTT-o l)exp[ ] 



I6bf 



32^2 



, m + kf 9{bi + k)\ ^ 
+ ( ,.,9 l)exp[ — ^] 



\6bj 



32bf 



and 

b) using full isotropy 

in position space: 



in Fourier space: 



(A.4) 



■exp[--bfr'-] 
27 V2^ 9 (A.5) 

((27 - I6b]r^) cos[bir] - 9^;r sin[^ir]) , 



e{k) = eoi-y = qk-'- 
ko 

The variance from the above power law is then defined as: 

ds - Si - qkj' 

The allowed statistical deviations of s are imprinted in a co- 
variance matrix Dij - {6si Ssj), which is automatically ob- 
tained from the inverse of the Hessian matrix calculated with 
the maximum a posteriori approach as described in Section 
13.31 The probabilities of the spectral amplitudes s of the 
power spectrum can be approximated as a Gaussian, where 
we can neglect the normalisation factor, because it is inde- 
pendent of q and r: 



P{ds\q, r) cc e 



kSs'D-'Ss 



. Only a part of the power spectrum is used to fit the power 
law. Djj is then the projection (cut) of the full inverse Hessian 

matrix Dij obtained from Eq. [16] 

Dij =< 6si6sj >= ij/Dijip, 
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with the projection operator - {0, 1, ...1, 0} selecting 
the used spectral bins for the fit. Using the Bayesian theorem 
from Eg.fTSl 

P{q, r\6s) oc P{6s\q, r), 

where we already applied a uniform prior distribution 
P{r, q) - const. P{q, r\6s) has to be maximised to achieve 
the desired slope r within a maximum-a-posteriori approach. 
To numerically find the optimum, the gradients and the 2x2 
Hessian matrix are necessary: 

^ = -q\n[ki]kT'Drl(sj-qk-/) 

^-^'Hr/is,-qkJ^) 

= ^(^ln[^,]^r'D7.>ln[^,]^7'- 
-{MkAfk-'-DjIiSj-qk-")) 

^'^"^ = 2q\n[ki]kJ'-Djlk-/- -\n[ki\kJ'Djlsj. 



An expression for the second unknown q can be found ana- 
lytically. If the initial value of r is given, the corresponding 
initial value q can be calculated easily. This is necessary be- 
cause the first guess must be close to the maximum to attain 
convergence. 
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